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ABSTRACT 


Results  of  an  analytical  study  to  theoretically  predict  the  loads  to  be 
experienced  by  a  Lockheed  Model  1649  Super  Constellation  during  a  con¬ 
trolled  crash  are  presented.  Acceleration  time  histories  in  directions 
normal  and  parallel  to  the  fuselage  cabin  floor  are  given  at  three  positions 
along  the  length  of  the  fuselage  for  impact  velocities  of  140,  160,  180,  and 
200  feet  per  second.  Results  of  investigations  of  the  effects  of  variations 
in  important  problem  parameters  are  also  presented. 

It  is  concluded  that  during  the  initial  impact  at  180  feet  per  second, 
peak  normal  accelerations  of  11,  0,  and  -3  times  that  of  gravity  (G’s)  will 
be  developed  0.03  seconds  after  impact  at  Body  Stations  180,  682,  and  1176. 
Maximum  normal  and  longitudinal  accelerations  during  the  6  degree  ramp 
crash  will  cccur  at  0.  24  seconds.  Maximum  normal  accelerations  at  Body 
Stations  180,  682,  and  1176  will  be  -17,  8,  and  35  G’s  respectively.  Max¬ 
imum  longitudinal  accelerations  will  be  4  G’s. 

It  is  further  concluded  that  the  nose  of  the  airplane  will  bend  upwards 
10  inches  relative  to  the  center  section  of  the  fuselage  0. 14  seconds  after 
impact.  This  deflection  will  probably  be  of  sufficient  magnitude  to  exceed 
the  ultimate  strength  of  the  fuselage  above  the  cabin  floor.  If  the  fuselage 
should  fail,  all  analytical  results  beyond  the  time  of  failure  will  be  question¬ 
able  to  a  degree  dependent  on  the  type  of  failure  that  occurs. 
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1. 0  INTRODUCTION 


This  report  presents  the  results  of  an  analytical  study  to  theoretically 
determine  the  loads  to  be  experienced  by  a  Lockheed  Model  1649  Super  Con¬ 
stellation  during  a  controlled  crash.  A  specific  requirement  of  this  study  was 
the  calculation  of  acceleration  time  histories  in  directions  parallel  and  normal 
to  the  cabin  floor  at  Fuselage  Stations  180,  682,  and  1176.  This  was  done  for 
each  of  fou*'  initial  impact  velocities:  140,  160,  180,  and  200  feet  per  second. 
Data  of  this  type  can  be  used  to  calculate  the  forces  that  would  be  experienced 
during  the  crash  by  passengers,  crew,  or  cargo  located  in  the  airplane. 

The  determination  of  the  loads  in  an  aircraft  during  a  marginally 
survivable  crash  is  a  complex  problem  and  depends  upon  many  factors.  Among 
the  more  important  of  these  factors  are  the  initial  velocity  of  the  airplane,  its 
attitude  relative  to  the  ground  and  to  its  flight  path,  its  geometry  and  weight, 
the  type  of  surface  on  which  it  impacts,  and  the  deformation  characteristics 
of  its  structure.  If  an  analytical  method  could  be  found  to  define,  with  a  reason¬ 
able  degree  of  accuracy,  marginally  survivable  crash  loads  for  different  types 
of  aircraft,  information  obtained  would  be  invaluable  in  the  development  of 
equipment  for  the  safety  of  passengers  and  crew.  Included  in  this  equipment 
would  be  such  items  as  airplane  structure,  passenger  seats,  safety  harnesses, 
shock  absorption  systems,  and  cargo  restraint  systems.  The  study  reported 
in  this  document  represents  an  attempt  to  develop  such  an  analytical  method. 

The  airplane  of  this  analysis  is  a  Lockheed  Model  1649  Super  Constel¬ 
lation  (Fig.  3).  Its  configuration  is  identical  to  that  to  be  used  by  the  Federal 
Aviation  Agency  in  an  actual  crash  test  of  the  Lockheed  1649  currently  scheduled 
for  August  1964.  The  specific  conditions  of  this  test  are  as  follows  (Ref.  2): 

•  The  airplane  is  to  impact  nose  first  into  a  6  degree  dirt  ramp  at  a  velocity 
of  approximately  180  feet  per  second. 

•  At  impact,  the  longitudinal  axis  of  the  fuselage  is  to  be  horizontal  to  the 
ground  at  the  base  of  the  ramp  and  the  wings  are  to  be  level. 

•  The  yaw  angle  of  the  airplane  relative  to  its  flight  path  is  zero. 

•  The  wing  lift  is  completely  destroyed. 

•  The  landing  gear  and  wings  outboard  of  Wing  Station  484  are  to  have  com¬ 
pletely  separated  from  the  airplane. 

These  initial  test  conditions  are  considered  to  be  of  a  marginally  sur¬ 
vivable  nature,  and  simulate  a  crash  condition  such  as  might  be  experienced  in 
an  aborted  takeoff. 

It  is  to  be  emphasized  that  the  results  of  this  report  apply  only  to  the 
airplane  as  it  passes  over  the  6  degree  ramp.  At  a  distance  of  approximately 
200  feet  from  the  initial  impact  point  the  ramp  slope  increases  sharply.  At 
this  point  these  results  are  no  longer  valid.  It  is  also  be  emphasized  that  this 
analysis  assumes  the  portion  of  the  airplane  above  the  cabin  floor  remains 
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intact  during  the  6  degree  ramp  crash.  If  this  part  of  the  fuselage  should  fail 
at  any  time  during  the  crash,  the  analytical  results  beyond  that  time  would  be 
questionable.  How  much  the  accuracy  of  the  results  would  be  affected  would 
depend  largely  on  the  type  of  failure  that  occurred.  For  example,  if  the  fuse¬ 
lage  forward  of  the  leading  edge  of  the  wing  separated  from  the  airplane,  the 
physical  characteristics  of  the  airplane  would  change  substantially,  and  ana¬ 
lytical  data  for  times  after  this  failure  would  not  be  reliable.  This  is  covered 
in  additional  detail  in  Section  4.  2. 1. 

The  discussion  of  this  report  is  divided  into  five  sections.  The  first 
section  explains  the  basic  equations  used  in  the  analysis.  These  equations  are 
derived  in  three  steps.  The  airplane  is  first  treated  as  a  rigid  body  free  to 
pitch  and  to  translate  horizontally  and  vertically  in  space.  Next,  the  vertical 
bending  flexibility  of  the  fuselage  is  included,  and,  finally,  equations  for  the 
longitudinal  vibration  of  the  fuselage  are  obtained.  Also  included  in  the  first 
section  is  a  brief  discussion  of  the  digital  program  used  for  the  solution  of  the 
mathematical  equations. 

The  second  section  of  this  report  discusses  the  computation  of  the 
basic  problem  input  data.  The  third  section  presents  the  major  results.  In 
the  first  part  of  this  section,  acceleration  time  histories  are  given  for  both  the 
rigid  body  analyses  and  for  those  in  which  the  vertical  bending  and  longitudinal 
flexibilities  of  the  fuselage  are  considered.  Because  of  the  large  number  of 
variables  involved  in  the  analyses  and  the  difficulties  of  mathematically  repre¬ 
senting  or  experimentally  determing  these  variables,  a  parametric  study  is 
described  in  the  second  part  of  the  third  section. 

The  important  conclusions  of  this  theoretical  study  and  recommenda¬ 
tions  for  future  work  are  given  in  the  final  sections  of  this  report. 

As  mentioned  previously,  the  calculation  of  crash  loads  is  complex 
and  necessarily  involves  many  approximations.  For  this  reason,  if  the  correct 
trends  of  the  acceleration  time  histories  are  predicted,  and  if  the  peak  accel¬ 
eration  values  and  times  at  which  they  occur  do  not  deviate  more  than  25  per¬ 
cent  from  the  experimental  data,  it  is  believed  good  agreement  between  ana¬ 
lytical  and  experimental  results  will  have  been  realized.  As  experience  is 
gained,  it  is  expected  that  the  accuracy  of  the  analysis  could  be  improved. 
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2.  0  BASIC  THEORY 


This  section  describes  the  mathematical  equations  used  to  determine 
the  motion  of  the  airplane  during  the  crash.  It  is  divided  into  three  main  sub¬ 
sections.  The  first  subsection  discusses  the  rigid  body  equations;  that  is, 
those  equations  in  which  the  vertical  bending  and  longitudinal  flexibilities  of  the 
airplane  are  not  considered.  The  second  gives  those  equations  which  account 
for  the  effects  of  airplane  vertical  bending  ;ind  longitudinal  vibrations.  The 
third  subsection  briefly  describes  the  techi  ique  used  in  the  solution  of  the 
mathematical  equations. 

2. 1  THE  RIGID  BODY  EQUATIONS 

2.1.1  Assumptions 

In  the  derivation  of  the  rigid  body  equations,  the  following  assumptions 
were  made; 

•  The  crash  is  symmetric.  That  is,  the  airplane  is  free  to  pitch  and  trans¬ 
late  vertically  and  horizontally  during  and  after  initial  impact,  but  cannot 
roll,  yaw,  or  tra  slate  laterally. 

•  At  initial  impact  the  wings  are  level  to  the  surface  of  the  ramp  and  the  yaw 
angle  of  the  airplane  relative  to  its  forward  motion  is  zero. 

•  The  wing  lift  is  completely  destroyed. 

•  The  bottom  of  the  fuselage  is  made  up  of  a  finite  number  of  nonrestoring, 
nonlinear  springs  in  parallel  with  nonreversible  linear  dampers. 

2. 1. 2  The  Equations 

The  following  are  the  equations  used  to  describe  the  rigid  body  motion.  A 
drawing  of  the  coordinate  system  is  shown  in  Fig.  1,  and  the  definitions  of  sym¬ 
bols  are  given  in  the  list  of  symbols. 

The  x-y  coordinates  of  the  ith  point  on  the  airplane  are: 

Xj  -  x  +  l  ■  cos  a  +  hj  sin  a  (1) 

Vj  ~  y  -  l  j  sin  a  +  h-  txxs  a  (2) 

The  normal  distance  from  the  ith  point  on  the  airplane  to  the  ramp  is 
therefore: 


Vj  -  Xj  sin  >'*  j  +  Vj  cos  r‘  j  (3) 

Note  that  the  ramp  can  have  two  angles  (Fig.  1)  with  the  change  in  ramp  slope 
occurring  at  the  origin  of  the  coordinate  system.  Note  also  that  according  to 
the  chosen  sign  convention  the  ith  point  is  in  contact  with  the  ramp  only  if  v.  is 


FIG.  1  SIGN  CONVENTION  USED  IN  DERIVATION  OF  EQUATIONS 

greater  than  zero.  Vi,  then,  can  be  used  directly  to  determine  the  deflection 
and,  therefore,  the  force  in  the  nonrestoring  springs  and  dampers  representing 
the  bottom  of  the  fuselage. 

Each  structural  spring  along  the  base  of  the  fuselage  is  assumed  to  be 
in  series  with  a  spring  representing  the  ground.  The  total  deflection  of  these 
two  springs  at  the  ith  station  must  then  equal  vi  for  values  of  Vi  greater  than 
zero. 


v,  =  ♦  et.  (4) 

Because  the  springs  are  in  series  the  force  in  each  of  these  springs  must  be 
equal  at  all  times. 


(5) 


A  viscous  damper  to  represent  the  combined  damping  characteristics 
of  the  ground  and  structure  can  be  placed  in  parallel  with  each  ith  spring. 
Again,  for  values  of  v*  greater  than  zero,  the  damping  force  is  then: 


VC- 

1  t 


(6) 


The  force  acting  on  the  airplane  normal  to  the  ground  surface  at  the  ith  location 
is: 


(7) 
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The  total  forces  acting  on  the  airplane  parallel  to  the  ground  surface 
are  considered  to  consist  of  friction  and  plowing  forces.  The  friction  force  at 
the  ith  station  is: 


Vwfty  (8) 

The  plowing  force  is  assumed  to  he  a  constant  force  acting  at  any  time  the  air¬ 
plane  is  in  contact  with  the  ramp.  This  is  not  correctly  true  since  the  plowing 
area  of  the  airplane  is  related  to  the  deformation  of  the  ground  and  to  the  air¬ 
plane  under  structure.  However,  because  of  the  difficulties  of  analytically  de¬ 
termining  exact  ground  and  airplane  understructure  deformation,  it  is  believed 
the  assumption  of  a  constant  plowing  force  is  adequate.  This  is  discussed  in 
more  detail  in  Section  4.  3.  5. 


The  plowing  force  is  then: 
C„ 


(9) 


Calculating  the  components  in  the  x  and  y  directions  of  forces  normal  and 
parallel  to  the  ramp: 


Fx  ~  ~  !  [  fN{  sin  6 1 +  F/ ■cos  ®i :  ]~FPcosei 


(10) 


Fy-  -z"_  t  [  F/I .  cos  0  ■  -  Ff.  sin  6  ,•  J  +  Fp  sine  J 


(11) 


Note  that  the  ramp  angle  at  the  nose  is  used  with  the  plowing  force.  The  mo¬ 
ment  arms  of  F  and  F,  about  the  airplane  center  of  gravity  are: 


•i  *i 

Rxy=h-  cos  a  -l -sin  a 

Ry .  =  h-  sin  a  +  £  ,•  cos  a 


And  the  total  moment  about  the  airplane  center  of  gravity  is: 

-  Z  n  t  Rx  (pN  ■  sin  6  i  +  Ff  cos  6  i}  + 

a  j  -  i  L  i  i  i 

Ry.fFff  cose  j  -  Ff'Sine  j)  J 


-  Rv  Fncos  6,  -  Rv  Ff-sine, 
Xj  p  i  p  i 


Setting  F^,  F  ,  and  Ma equal  to  the  inertia  forces  in  the  system  gives: 


Mx  =  F. 


My  -  Fy  +  MG 


I'd 


(12) 

(13) 


(14) 


(15) 

(16) 
(17) 
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Using  the  above  equations  to  determine  x,  y,  and  a  ,  the  relative 
accelerations  parallel  and  normal  to  the  fuselage  floor  at  the  jth  point  can  be 
calculated: 

••  ••  - 

xj  =  xcos  a.  -  y  sina  -  t  j  a  2  +  hj  a  +  G  sin  a  (18) 

pj  -  —‘i  sin  a  —  y  cos  a  +  hj  a  2  +  lj  &+  G  cos  a  (19) 

The  first  two  terms  of  each  of  the  above  equations  are  the  components 
of  the  x  and  y  accelerations  acting  parallel  and  normal  to  the  fuselage  floor. 
The  third  terms  account  for  centrifugal  force  effects,  and  the  fourth  for  accel¬ 
erations  caused  by  pitch  of  the  airplane.  The  final  terms  remove  the  effect 
of  gravity  making  the  accelerations  relative  rather  than  absolute. 

Equations  (1)  through  (19),  then,  comprise  the  basic  rigid  body  equa¬ 
tions  for  the  symmetrical  crash  of  an  airplane  into  a  wedge-shaped  ramp. 

2.  2  THE  FLEXIBLE  BODY  EQUATIONS 

In  the  derivation  of  the  flexible  body  equations,  the  assumptions  of 
Section  2. 1. 1  were  used.  As  the  crash  was  considered  to  be  symmetric,  the 
only  airplane  flexibilities  considered  were  those  in  the  vertical  bending  and 
longitudinal  directions.  These  will  be  discussed  individually. 

2.  2. 1  Vertical  Bending  of  the  Fuselage 

The  displacement  of  the  fuselage  in  its  first  vertical  bending  mode  is 
represented  by  the  single  generalized  coordinate  “q.  ”  The  displacement  nor¬ 
mal  to  the  fuselage  longitudinal  axis  of  the  ith  point  on  the  fuselage  is  then  4>  jq 
where  4»j  is  the  mode  shape  weighting  function  as  given  in  Fig.  6. 

Calculating  the  kinetic  and  potential  energies  of  the  airplane  in  its 
first  vertical  bending  mode,  and  differentiating  according  to  the  Lagrangean 
equation 


gives  the  equation  of  motion  for  the  vertical  bending  degree  of  freedom 

=  Fq  (21) 

where  E  is  the  generalized  mass  of  the  fuselage  in  its  first  vertical  bending 
mode  and  is  defined  by  the  equation 


The  right  side  of  Equation  (21)  consists  of  two  terms.  The  first  is 
the  generalized  stiffness  of  the  vertical  bending  degree  of  freedom.  The  sec- 
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ond  term  is  the  generalized  force  and  is  determined  by  the  equation 

n 

V1  1=1  °iCOS(  ’i"  Q  )+F/;  *isin(  Qi“  aJ]  0*3) 

+  Fp  C  j  sin  (  9  j  -  a  ) 

v.  must  now  be  redefined.  Taking  q  as  positive  nose  down 


v,  '  x;  s/n  e  i t  >7  cos  0  ,  +  cos  r?,-  a  )- i-  q  (24) 

The  accelerations  normal  and  parallel  to  the  cabin  floor  at  the  jth  point  are: 


% 

1(j  -  x  cos  a  —  y  sin  a  — 

+  hi  a  +  G  sin  a 

1 

(25) 

y.  =  _x  sin  a  -  y  cos  a  +  hj 

a  ^  +  l  ■  *a  +  G  cos  a  -  £  ■  q 

1  J 

(26) 

With  the  exception  ot  Equation  (3),  Equations  (1)  through  (17)  from 
Section  2. 1.2  together  with  (21)  through  (26)  above  are  the  equations  describing 
the  symmetrical  crash  of  an  airplane  having  fuselage  vertical  bending  flexibility. 

2.  2.  2  Longitudinal  Fuselage  Vibrations 

In  addition  to  the  assumptions  of  Section  2. 1. 1,  the  following  assump¬ 
tions  were  used  in  the  derivation  of  the  equations  for  fuselage  longitudinal  vi¬ 
brations. 


1.  The  fuselage  is  a  continuous  beam  with  constant  weight,  density, 
and  bulk  modulus  of  elasticity  along  its  length. 

2.  The  continuous  beam  representing  the  fuselage  is  in  a  fixed-free 
end  condition.  That  is,  the  nose  of  the  airplane  is  not  free  to  displace  longi¬ 
tudinally  relative  to  its  rigid  body  displacement. 

3.  The  longitudinal  vibrations  are  dependent  only  on  the  longitudinal 
accelerations  experienced  by  the  airplane  during  the  crash  and  are  not  affected 
by  forces  or  accelerations  normal  to  the  longitudinal  axis  of  the  fuselage. 

4.  The  total  longitudinal  acceleration  of  the  airplane  is  the  sum  of 
its  rigid  and  flexible  body  accelerations. 

These  assumptions  are  severe,  and  were  made  to  allow  approximate 
answers  to  be  obtained  in  a  reasonable  period  of  time.  A  more  accurate  mathe¬ 
matical  representation  could  be  derived,  but  the  complexity  of  the  equations 
would  be  formidable.  Since  results  have  shown  that  the  effects  of  longitudinal 
vibrations  of  the  fuselage  are  small  relative  to  the  rigid  body  accelerations, 
and  would  probably  be  even  smaller  if  more  exact  mathematical  equations  were 
used,  the  theory  of  this  report  is  believed  to  be  adequate. 
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C  D 


FIG.  2  FORCES  ACTING  ON  ELEMENT  OF  BAR  IN  LONGITUDINAL  MOTION 


Assuming  a  beam  of  constant  weight,  density,  and  bulk  modulus,  the 
forces  acting  on  an  element  of  the  beam  can  be  written.  Consider  a  small  ele¬ 
ment  CD  of  length  Ax  (Fig.  2)  and  let  the  cross-sectional  area  of  the  beam  be 
A.  The  unit  stress  on  the  face  through  C  is  a  x  and  the  stress  on  the  face 
through  D  is  therefore  ox  +  (3ax/3*)Ax.  Letting  the  longitudinal  displace¬ 
ment  of  the  element  be  u,  we  have  from  Newton’s  second  law  of  motion: 
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Equation  (30)  then  becomes: 


x 


.s 


2 


The  following  are  the  boundary  conditions: 


F(t) 


And 


~s  (L,  i ) 


(32) 


where  L  is  the  x  coordinate  of  the  beam  at  its  free  end  and  F(t)  is  the  longi¬ 
tudinal  force  caused  by  the  rigid  body  motions. 

Equation  (32)  can  be  solved  by  a  numerical  method.  This  is  explained 
more  completely  in  Ref.  3.  The  general  procedure  is  to  define  the  impact  force 
function  F(t)  at  the  nose  of  the  airplane  (x  =  0).  The  initial  "tress  F(t)/A  is 
considered  to  travel  along  the  bar  as  a  compressive  wave.  Upon  reaching  the 
aft  or  free  end  of  the  fuselage  it  reflects  back  as  a  tensile  stress  wave.  When 
this  reflected  wave  reaches  the  nose  or  impact  end  of  the  fuselage  it  again  re¬ 
flects  back  through  the  fuselage  as  a  compressive  stress  wave  with  reduced 
intensity.  Looking  at  a  particular  point  xj  along  the  length  of  the  fuselage  it 
initially  takes  time 

*o ~  (33) 


for  the  stress  wave  to  reach  points  xj.  It  also  takes  time 
2L  -  x 

t  _ L 

1  ~  (34) 


for  the  stress  wave  to  reach  point  xj  after  reflecting  back  from  the  free  end 
x  =  L.  This  process  is  repeated  with  period  2L/s,  and  the  magnitude  de¬ 
creases  each  time  by  a  factor  r.  The  stress  at  xj,  then,  can  lx*  obtained  by 
algebraically  summing  the  original  and  reflected  stress  waves.  The  mathe¬ 
matical  expression  of  the  resultant  stress  can  lx?  written: 


’.V  (*j-  ” 


1 

~A 


/  r 


Fii-Ia)  -  v  F(t-t 


m  0 


B 


(35) 


where: 


• 

<A  lo  *  "r 

(30) 

1 B  U  "  m 

(37) 

• 

and  B  and  ■  are  step  functions  such  that: 

0  For  0  -  r  - 

1.0  For  t  ■  tA 

(3«) 

(39) 


>  _  o  F or  0  -  l  -iB 

‘  IS)  Fur  t  >  tB 

The  first  term  in  (35)  within  the  range  of  the  summation  sign  repre¬ 
sents  the  compressive  waves,  and  the  second  term  is  the  contribution  to  the 
total  stress  of  the  tensile  waves. 

Having  obtained  the  stress  function  x  (Xj ,  t),  the  acceleration  at 
station  x-  can  be  calculated  by  differentiating  x  with  respect  to  x  at  the  jth 
point  anti  substituting  in  Equation  (27). 

°i  =—  (  -7-)  i  (40) 

The  total  longitudinal  acceleration  at  any  point  along  the  fuselage  is 

then: 

(ii) 

by  Assumption  (4) 

2.  3  SOLUTION  OF  THE  EQUATIONS 

The  differential  equations  of  motion  derived  in  Sections  2. 1  and  2.  2  are 
of  second  order  and  are  nonlinear.  For  this  reason,  the  technique  of  numerical 
integration  was  used  to  solve  them.  The  solution  program  employed  was  the 
dynamic  analyzer  (DYANAj  program  described  in  Ref.  7,  This  program  utilizes 
the  Runge-Kutta  method  of  numerical  integration,  a  description  of  which  can  be 
found  in  Refs.  11  and  13.  Essentially,  in  numerical  integration,  the  accelera¬ 
tion,  velocity,  and  displacement  of  the  coordinates  of  the  system  are  computed 
over  a  given  time  interval.  These  computed  qualities  are  then  used  as  initial 
conditions  for  the  next  time  interval  in  the  evaluation  of  acceleration,  velocity, 
and  displacement  at  the  end  of  that  time  interval.  Generally  speaking,  for  a 
given  time  range  decreasing  the  time  interval  of  integration  increases  the 
accuracy  of  the  answer,  hi  the  solution  of  the  equations  of  Sections  2.  1  and 
2.  2,  a  time  range  of  0.  9  seconds  was  selected,  and  integrations  were  performed 
at  intervals  of  0.  001  seconds  within  that  range.  Answers  were  pointed  at  incre¬ 
ments  of  0.  02  seconds. 
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3.0  BASIC  INPUT  DATA 

This  section  discusses  the  calculation  of  the  basic  problem  input  data. 

A  drawing  of  the  Model  1649  is  given  in  Fig.  3. 

3.  1  AIR  PLANK  WEIGHT  AND  INERTIA 

The  airplane  weight  properties  were  estimated  from  data  supplied  in 
both  Lockheed  and  the  Federal  Aviation  Agency.  Fig.  1  shows  the  weight  dis¬ 
tribution  used  in  the  analysis.  The  separation  from  the  airplane  of  the  wings 
outboard  of  Wing  Station  484  has  been  accounted  for  in  this  listing. 

3.  2  AIRPLANE  STRUCTURAL  STIFFNESS  DATA 

3.2.  1  Flexibility  of  the  Fuselage  Understructure 

During  impact  of  the  airplane  into  the  ramp,  energy  is  absorbed  by 
the  crushing  of  the  understructure  of  the  fuselage.  It  is  important,  then,  that 
in  a  crash  loads  analysis  that  the  crushing  characteristics  of  the  fuselage  under¬ 
structure  be  represented. 

The  bottom  of  the  fuselage  is,  of  course,  a  continuous  stmeture.  It 
was  necessary  for  caleulative  purposes,  however,  to  consider  it  as  consisting 
of  a  finite  number  of  springs  in  parallel  along  the  length  of  the  fuselage.  A 
total  of  20  were  used,  and  a  listing  of  their  properties  is  given  in  Fig.  3.  With 
reference  to  the  sample  stiffness-deflection  curve  shown  at  the  top  of  Fig.  3. 
the  springs  were  considered  to  be  elastic  (restoring!  within  the  range  of  j . 
When  the  deflection  of  the  spring  increased  Itcvond  j.  however,  the  structure 
was  considered  to  be  permanently  deformed  (plastic)  and  the  spring  was  theie- 
lorc  not  restoring.  The  downward  slope  of  the  stiffness  curve  between  i  and 
2  represents  the  decrease  in  effective  spring  stiffness  because  of  exceedence 
of  the  elastic  limit.  The  upward  slope  of  the  curve  beyond  _>  accounts  for  the 
increase  because  of  the  stiffen  structure  of  the  cabin  floor  and  wing. 

Basic  information  on  1(549  fuselage  frame  and  bulkhead  spacing  and 
construction  was  furnished  by  Lockheed.  This  was  used  to  determine  the  spring 
properties  of  Fig.  3.  The  general  procedure  was  io  div  ide  the  fuselage  into 
sections  of  varying  length  based  on  geometry  and  location  of  bulkheads.  The 
deflection  per  unit  load  of  the  frames  and  bulkheads  within  each  section  was 
then  calculated  using  the  methods  of  Ref.  12.  This  gave  the  stiffness  of  the 
section  up  to  that  deflection  at  which  the  frame  or  bulkhead  yields  due  to  elas¬ 
tic  instability.  This  stiffness  is  shown  as  Kj  in  the'  drawing  ol  Fig.  3.  The 
stiffness  of  the  sections  for  deflections  beyond  the  elastic  limit  was  estimated. 

1!  the  geonietiw  of  each  section  did  not  change  as  the  understructurc  was 
crushed,  it  would  be  expected  that  once  the  elastic  limit  was  exceeded,  the 
crushing  force  would  remain  relatively  constant  as  the  deflection  increased. 
However.  sinc<  .  as  the  understructurc  of  the  1619  crushes,  additional  mate¬ 
rial  picks  up  the  load,  it  is  believed  the  estimated  stiffness  values  give  a  good 
approximation  of  the  force-deflection  characteristics  of  the  stmeture. 

In  the  calculation  of  this  stiffness  data,  the  contribution  of  the  nacelles 
was  not  considered  to  be  significant,  'ibis  may  not  be  strictly  true  as  the  oil 


FIG.  3  THE  LOCKHEED  MODEL  1649  SUPER  CONSTELLATION 
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cooler  air  intake  duets  on  the  bottom  of  the  inboard  nacelles  extend  about  one 
foot  below  the  fuselage.  It  is  believed,  however,  that  the  nacelle  structure  is 
relatively  weak,  and  will  be  separated  from  the  airplane  when  it  contacts  the 
propeller  barriers  prior  to  ramp  impact  (Ref.  2). 

d.  2.  2  Body  Flexibility 

The  vertical  bending  stiffness  (El)  of  the  fuselage  for  the  calculation 
of  the  fuselage  vertical  tending  mode  was  furnished  by  Loekheed.  A  listing  of 
this  data  is  given  in  Fig.  4. 

B.S.  B.S. 


PANEL 

i  n 

I 

PANEL 

NO. 

BODY 

STATION 

WEIGHT 

LBS 

E.l.  x  10  ,u 

lb  ir 

1 

NO. 

BODY 

STATION 

WEIGHT 

LBS 

TXTio 10 

LB  IN^ 

1 

138.4 

514 

1.0 

14 

783.1 

2331 

80.0 

2 

193.4 

2471 

3.0 

15 

83b. S 

1852 

79.0 

3 

239.5 

2359 

10.0 

16 

884.1 

16S2 

63.0 

4 

284.5 

i872 

17.5 

17 

933.0 

1656 

58.0 

5 

334.0 

1844 

25.0 

15 

984.1 

1685 

54.0 

6 

384.6 

1809 

33.0 

15 

1034.6 

1579 

40.0 

i 

435.6 

I  707 

I I  Lj 

38.0 

20 

1083.5 

1542 

33.0 

8 

433./ 

4815 

39.5 

21 

1130.4 

503 

26  u 

9 

534.4 

4805 

47.0 

22 

1182.9 

377 

18.0 

10 

584.7 

i8S7 

5S.5 

23 

1237.1 

343 

12.5 

11 

616.7 

2903 

75.0 

24 

1269.8 

o-'n 

i  Jj 

10.0 

V  .,0 

648  0 

75748 

RIGID 

25 

1330.9 

255 

7.5 

12 

681.7 

2263 

RIGID 

26 

1391.9 

3254 

5.0 

13 

731  4 

1940 

RIGID 

27 

1428  3 

225 

1.0 

TOTAL  582.0  124  500 


FI G.  4  BASIC  WEIGHT  AND  STIFFNESS  DATA  -  THE  LOCKHEED  MODEL  1649 


:i.:t  fuselage  v  f  r  ti  c  a  I  bending  mode 

l'o  obtain  the  frequency  and  nuxle  shape  of  the  fuselage  tirst  vertical 
bending  nuxle,  the  fuselage  was  divided  into  three  principal  sections:  that  sec¬ 
tion  of  the  fuselage  forward  of  and  cantilevered  from  the  wing  front  spar,  the 
lusealge  aft  of  and  cantilevered  from  the  wing  rear  spar,  and  the  section  of  the 
luselage  between  the  wing  spars.  The  fore  and  aft  fuselage  sections  were  con¬ 
sidered  to  be  flexible  while  the  center  body  was  treated  as  a  rigid  member. 

The  first  cantilevered  vertical  bending  mode  of  the  fore  and  aft  fuselage  were 
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determined  using  the  influence  coefficient  method  described  in  Chapter  4  of 
Ref.  14.  These  vertical  bending  modes  of  the  fore  and  aft  bodies  were  then 
coupled  with  rigid  body  pitch  and  vertical  translation  by  means  of  Lagrange’s 
Equation  (20).  The  Rayleigh-Ritz  method  Jso  described  in  Ref.  14  was  then 
used  to  compute  the  first  unrestrained  vertical  bending  mode  of  the  fuselage. 
The  normalized  mode  shape  is  given  in  Fig.  0. 

3.  4  OTHER  DATA 

Information  on  the  calculation  of  other  parameters  such  as  fuselage 
plowing  area,  friction  coefficient,  ground  stiffness,  etc. ,  is  contained  in  Sec¬ 
tion  4.  3. 


1 

BODY 

STA. 

FUSELAGE 

SECTION 

LENGTH 

FUSELAGE 

DIAMETER 

SPRING 

LOCATION 

K1 

K 

2 

6 

1 

S2 

TANGENT 

*1 

□ 

r 

INCHES 

INCHES 

B.S. 

LB./IN. 

LB./IN./IN. 

!  i 

100 

54 

130 

12000 

5000 

2 

8 

-1165 

1165 

■I;  j® ;  1 

80 

100 

230 

9600 

4000 

3 

20 

-  306 

306 

®i;  |®  \  1 

60 

120 

290 

6600 

3000 

27 

-  150 

15l 

4 

320-  380 

130 

350 

31 

-  129 

129 

5 

380-  440 

138 

410 

35 

-  112 

1 

2 

6 

440-  500 

140 

470 

1 

500-  560 

530 

8 

560-  620 

6 

0 

590 

6600 

3000 

-  1 

2 

1 

2 

9 

620-  682 

6 

0 

c 

652 

12400 

5600 

-  2 

2 

1250 

10 

682-  744 

62 

714 

12400 

5600 

-  2 

2 

1250 

11 

744-  605 

61 

775 

12200 

5600 

-  2 

12 

1250 

12 

805-  862 

57 

832 

6850 

2900 

-  125 

125 

13 

862-  919 

57 

14 

0 

889 

6850 

2900 

-  125 

12c 

14 

919-  989 

70 

13 

0 

959 

7000 

3500 

-  110 

no 

15 

989-1059 

120 

1029 

35 

-  110 

no 

16 

1059-1129 

114 

1099 

29 

-  135 

135 

1/ 

1129-.  199 

106 

1169 

22 

-  184 

184 

18 

1199-1269 

90 

1239 

3 

15 

-  290 

290 

19 

1269-1339 

7 

D 

74 

1309 

70 

00 

3& 

30 

2 

9 

-  500 

500 

20 

1339-1410 

71 

66 

1379 

35 

00 

1750 

_ 

2 

4 

-  875 

875 

FIG.  5  FUSELAGE  UNDERSTRUCTURE  STIFFNESS  DATA 
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4. 0  RESULTS 


This  section  is  divided  into  three  main  subsections.  The  first  presents 
the  results  of  the  rigid  body  analyses,  and  the  second  those  of  analyses  consid¬ 
ering  the  longitudinal  and  vertical  bending  flexibilities  of  the  airplane.  The 
third  subsection  discusses  the  effects  of  variations  in  important  analysis  para¬ 
meters. 

4.  i  RIGID  BODY  ANALYSES 

Because  it  was  believed  from  the  beginning  of  these  analyses  that  the 
rigid  body  motions  would  dominate,  it  was  decided  to  present  the  rigid  body  re¬ 
sults  separately  from  those  results  where  the  longitudinal  and  vertical  bending 
flexibilities  of  the  fuselage  were  included.  In  the  rigid  body  analysis  the  only 
motions  considered  were  those  of  airplane  center  of  gravity  horizontal  and 
vertical  translation,  and  airplane  pitch  about  the  center  of  gravity. 

Figs.  7  through  10  show  the  acceleration  versus  time  histories  during 
the  6  degree  ramp  crash  determined  using  the  rigid  body  equations  of  Section 
2. 1.  Accelerations  normal  and  parallel  to  the  fuselage  cabin  floor  at  Body 
Stations  180,  682,  and  1176  are  given  for  initial  impact  velocities  of  140,  160, 
180,  and  200  feet  per  second.  A  summary  plot  showing  peak  normal  accelera¬ 
tions  for  each  velocity  and  time  at  which  the  peaks  occur  is  contained  in  Fig. 

11. 


The  basic  mechanism  of  the  crash  is  illustrated  in  Fig.  12.  In  this 
figure,  airplane  position  is  plotted  against  normal  acceleration  at  the  center 
of  gravity  and  airplane  velocity  parallel  to  the  ramp  surface.  At  time  zero  the 
airplane  impacts  the  ramp  approximately  30  feet  from  the  ramp’s  base  at  180 
feet  per  second.  The  longitudinal  axis  of  the  airplane  is  parallel  to  the  horiz¬ 
ontal.  At  impact  a  force  is  quickly  developed  at  the  nose  of  the  fuselage  which 
causes  the  airplane  to  pitch  up  sharply,  and  this  is  reflected  in  the  initial  ac¬ 
celeration  peaks  at  a  time  after  impact  of  0.  08  seconds.  The  aircraft  fuselage 
is  now  pitching  into  the  ramp,  and  finally  at  a  time  of  approximately  0.24 
seconds,  the  entire  fuselage  forward  of  the  trailing  edge  of  the  wing  impacts. 

At  this  point,  the  peak  accelerations  of  the  6  degree  ramp  crash  occur.  Be¬ 
cause  of  the  upward  curvature  of  the  aft  portion  of  the  fuselage,  the  aft  body 
does  not  contact  the  ramp  surface  until  after  the  peak  accelerations  have  passed. 
The  crushing  of  the  aft  fuselage  is  reflected  by  the  gradual  decrease  of  the  ac¬ 
celeration  from  its  maximum  value. 

Note  that  the  velocity  of  the  airplane  after  0.  7  seconds  is  only  25  feet 
per  second  less  than  at  impact.  It  is  estimated  that  it  will  take  approximately 
0.7  seconds  after  initial  impact  for  the  airplane  nose  to  reach  the  top  of  the 
6  degree  ramp. 

As  can  be  seen  by  a  study  of  the  acceleration  plots,  the  behavior  of 
the  airplane  after  impact  is  much  the  same  for  all  velocities.  The  most  signifi¬ 
cant  differences  appear  to  be  that  as  the  velocity  increases,  the  peak  accelera¬ 
tions  increase  and  the  time  at  which  they  occur  decreases.  This  trend  is  il¬ 
lustrated  clearly  by  Fig.  II. 
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INITIAL  CONDITIONS  -  RIGID  BODY 


SPEED 

140  FT. /SEC. 

C.G.  IiOCATION-  BODY  STA. 

682 

ANGLE  OF  IMPACT 

6  DEGREES 

PLOWING  COEFFICIENT 

80 LB.  /IN.2 

WEIGHT 

124,500  LBS. 

PLOWING  AREA 

700  IN.2 

INERTIA  ABOUT  C.G.  4. 64x10  JLB. /IN.  z 

FRICTION  COEFFICIENT 

.3 

NORMAL  ACCELERATION  -  G  UNITS 


FIG.  7  ACCELERATION  TIME  HISTORY  -  RIGID  BODY  -  140  FT  SEC  (continued) 
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LONGITUDINAL  ACCELERATION  -  G  UNITS 


SPEED 

100  FT.  /SEC. 

C.G.  LOCATION  -  BODY  STA. 

682 

ANGLE  OF  IMPACT 

G  DEGREES 

PLOWING  COEFFICIENT 

S0LB./IN.2 

WEIGHT 

124,500  LBS. 

PLOWING  AREA 

700  IN.2 

INERTIA  ABOUT  C.G. 

4.64xl09LB./IN.2 

FRICTION  COEFFICIENT 

.3 

FIG.  8  ACCELERATION  TIME  HISTORY  -  RIGID  BODY  -  160  FT/SEC 
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NORMAL  ACCELERATION  -  G  UNITS 


FIG.  8  ACCELERATION  TIME  HISTORY  -  RIGID  BODY  160  FT  SEC  (co  tinued) 
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LONGITUDINAL  ACCELERATION  -  G  UNITS 


INITIAL  CONDITIONS  -  RIGID  BODY 


SPEED 

ISO  FT. /SEC. 

C.G.  liOCATION  -  BODY  STA.  6H2 

ANGLE  OF  IMPACT 

‘i  DEGREES 

PIXIWING  COEFFICIENT  80LB./IN.2 

WEIGHT 

124,500  LBS. 

PI/HYING  AREA  700  IN.2 

INERTIA  ABOUT  C.G. 

4.!Mxl09LB./IN.2 

FRICTION  COEFFICIENT  .3 

FIG.  9  ACCELERATION  TIME  HISTORY  -  RIGID  BODY  -  180  FT /SEC 


A 
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NORMAL  ACCELERATION  -  G  UNITS 


FIG.  9  ACCELERATION  TIME  HISTORY  -  RIGID  BODY  -  180  FT/SEC  (continued) 
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li  .1  tV  all 

TIME  AFTER  IMPACT  -  SECONDS 

FIG.  JO  ACCELERATION  TIME  HISTORY  -  RIGID  BODY  -  200  FT/SEC 
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STA.  1176 


INITIAL  VELOCITY  -  FT/SEC. 


An  additional  interpretation  of  the  results  can  be  made  in  terms  of 
what  passengers  or  crew  seated  at  Body  Stations  180,  682,  and  1176  would  feel 
during  the  6  degree  ramp  crash.  For  example,  Figs.  7  through  10  indicate 
that  at  Body  Station  180,  a  passenger  at  initial  impact  would  be  thrown  forward 
and  downward  into  his  seat,  while  a  passenger  at  Body  Station  1176  would  be 
accelerated  forward  and  upward.  At  the  airplane  center  of  gravity,  Body 
Station  682,  the  passengers  would  experience  a  slight  downward  and  forward 
acceleration. 

At  peak  G’s  at  an  approximate  time  of  0.  22  seconds,  passengers  at  all 
stations  would  experience  maximum  forward  accelerations.  Those  passengers 
at  Body  Stations  682  and  1176  would  be  thrown  downward  into  their  seats  while 
those  in  the  nose  at  Body  Station  180  would  be  accelerated  upward  with  con¬ 
siderable  force.  It  would  appear  that  because  the  airplane  is  pitching  about  its 
center  of  gravity,  passengers  seated  at  or  just  forward  of  the  airplane  center 
of  gravity  would  experience  the  least  overall  forces  during  the  crash. 

4.2  FLEXIBLE  BODY  ANALYSES 

In  a  symmetric  crash,  as  is  being  considered  here,  fuselage  longi- 
tundinal  and  vertical  bending  vibrations  might  modify  the  rigid  body  results 
significantly.  Accordingly,  the  equations  of  Section  2.0  were  used  to  obtain 
acceleration  time  histories  similar  to  those  for  the  rigid  body  analysis.  In 
Figs.  13  through  16  accelerations  normal  and  parallel  to  the  fuselage  cabin 
floor  at  Body  Stations  180,  682,  and  1176  are  given  for  initial  impact  velocities 
of  140,  160,  180,  and  200  feet  per  second.  A  summary  plot  of  peak  normal  G’s 
at  the  airplane  center  of  gravity  versus  time  after  impact  is  shown  in  Fig.  17. 

In  general,  introduction  of  the  fuselage  vertical  bending  and  longi¬ 
tudinal  flexibility  into  the  problem  had  little  direct  effect  on  the  answers.  The 
rigid  body  accelerations  tended  to  control  the  results.  One  significant  difference 
that  is  noted  is  in  the  normal  acceleration  peaks  at  the  nose  and  tail  at  a  time 
after  impact  of  approximately  0.08  seconds.  Analysis  indicates  that  the  nose 
and  tail  at  initial  impact  tend  to  bend  upward  relative  to  the  relatively  rigid 
center  body  at  the  wing  fuselage  juncture.  This  does  not  significantly  affect 
the  normal  accelerations  at  the  center  of  gravity,  but  it  does  decrease  the 
initial  normal  accelerations  at  the  tail  and  causes  the  normal  acceleration  peak 
at  Body  Station  180,  which  occurs  0.  08  seconds  after  impact  in  the  rigid  body 
analysis,  to  shift  to  0.03  seconds. 

4.2.1  Fuselage  Failure 

Fig.  18  shows,  plotted  against  impact  velocity,  maximum  fuselage 
bending  deflections  at  the  nose  during  initial  impact.  All  maximum  deflections 
occurred  0. 14  seconds  after  impact.  An  estimate  of  the  ultimate  strength  of 
the  fuselage  indicates  that  a  bending  deflection  cf  approximately  5-8  inches  at 
the  nose  would  be  sufficient  to  fail  the  fuselage  in  the  neighborhood  of  the  for¬ 
ward  spar  of  the  wing  at  the  wing-body  juncture.  Therefore,  although  margin¬ 
al  at  the  lower  velocities,  it  would  appear  that  the  1649  fuselage  will  fail  during 
the  test.  Since  the  analysis  assumes  the  portion  of  the  fuselage  above  the 
cabin  floor  remains  intact  during  the  6  degree  ramp  crash,  analytical  results 
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beyond  the  time  of  failure  will  be  questionable  since  the  physical  characteristics 
of  the  airplane  will  change.  The  degree  to  which  the  accuracy  of  the  results 
is  affected  will  depend  on  the  type  of  failure  that  occurs. 

4.3  PARAMETRIC  VARIATIONS 

The  loads  experienced  by  an  aircraft  during  a  crash  are  dependent  on 
a  number  of  parameters,  many  of  which  cannot  always  be  determined  or 
analytically  represented  with  desired  engineering  accuracy.  The  most  imjx)r- 
tant  of  these  are  the  following: 

•  Impact  velocity 

•  Airplane  geometry,  weight,  and  inertia 

•  Ramp  shape 

•  Friction  force 

•  Plowing  force 

•  Damping  forces 

•  Ground  stiffness 

The  relative  effects  of  each  of  these  parameters  will  be  discussed 
individually  in  this  section.  All  analyses  accomplished  for  this  section  utilized 
only  the  rigid  body  equations. 

4. 3. 1  Impact  Velocity 

The  influence  of  initial  impact  velocity  has  already  been  discussed  in 
Section  4.  2.  Figs.  11  and  17  most  clearly  illustrate  how  peak  accelerations 
and  the  time  at  which  they  occur  vary  with  increasing  impact  velocity.  Gener¬ 
ally,  increasing  the  impact  velocity  causes  higher  peak  accelerations  at  an 
earlier  time  after  impact. 

4.3.2  Airplane  Geometry,  Weight,  and  Inertia 

Studies  of  the  effect  of  variations  of  airplane  weight  and  inertia  on 
crash  loads  were  not  accomplished  for  the  Lockheed  1649.  However,  some 
idea  of  the  magnitude  of  the  effects  of  airplane  geometry,  weight,  and  inertia 
can  be  gained  by  a  study  of  Fig.  19.  This  presents  the  results  of  an  analytical 
crash  loads  investigation  of  a  Curtiss  C-46,  an  airplane  with  a  gross  weight  of 
approximately  30  percent  that  of  the  Constellation.  *  The  initial  conditions  for 
the  C-4G  study  were  exactly  the  same  as  those  for  the  1649  investigation;  that 


This  C— 46  analysts  uas  act'omplished  at  an  earlier  date,  and  duplicates  an  experimental  test 
conducted  the  .V/lC/l  as  reported  in  Ref,  9. 
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LONGITUDINAL  ACCELERATION 


INITIAL  CONDITIONS  -  FLEXIBLE  BODY 
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FIG.  13  ACCELERATION  TIME  HISTORY  -  FLEXIBLE  BODY  -  140  FT/SEC 
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NORMAL  ACCELERATION  -  G  UNITS 


TIME  AFTER  IMPACT  -  SECONDS 

FIG.  13  ACCELERATION  TIME  HISTORY  -  FLEXIBLE  BODY  -  140  FT  SEC  (continued) 
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SPEED 

160  FT.  /SEC. 

C.G.  LOCATION- BODY  STA.  682  : 

ANGLE  OF  IMPACT 

6  DEGREES 

PLOWING  COEFFICIENT  80LB./IN.2 

WEIGHT 

124,500  LBS. 

PLOWING  AREA  700  IN.2 

INERTIA  ABOUT  C.G.  4. 64xl09LB.  /IN.  2 

FRICTION  COEFFICIENT  .3 

) 


FIG.  14  ACCELERATION  TIME  HISTORY  -  FLEXIBLE  BODY  -  160  FT/SEC 
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NORMAL  ACCELERATION  -  G  UNITS 


20 


: 
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TIME  AFTER  IMPACT  -  SECONOS 

FIG.  14  ACCELERATION  TIME  HISTORY  -  FLEXIBLE  BODY  -  160  FT  SEC  ( continued ) 


33 


LONGITUDINAL  ACCELERATION 


.2  .4  .6  .8 

TIME  AFTER  IMPACT  -  SECONDS 


FIG.  15  ACCELERATION  TIME  HISTORY  -  FLEXIBLE  BODY  -  180  FT  SEC 
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NORMAL  ACCELERATION  -  G  UNITS 


TIME  AFTER  IMPACT  -  SECONDS 

FIG.  15  ACCELERATION  TIME  HISTORY  -  FLEXIBLE  BODY  -  1 80  FT/SEC  (continued) 
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FIG.  16  ACCELERATION  TIME  HISTORY  -  FLEXIBLE  BODY  -  200  FT  SEC 


NORMAL  ACCELtRATION  -  G  UNITS 


TIME  AFTER  IMPACT  -  SECONDS 


FIG.  16  ACCELERATION  TIME  HISTORY  -  FLEXIBLE  BODY  -  200  FT/SEC  (confnutd) 
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FIG.  17  ACCELERATION  TIME  -  VELOCITY  SUMMARY  -  FLEXIBLE  BODY 
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INITIAL  VELOCITY  -  FT/SEC. 
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FIG.  18  MAXIMUM  BENDING  DEFLECTIONS  AT  NOSE 
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INITIAL  CONDITIONS 


SPEED  160  FT. /SEC. 

ANGLE  OF  IMPACT  6  DEGREES 

WEIGHT  37, *00  LBS. 

INERTIA  ABOUT  C.G.  1.5xl09LB. /IN. 2 
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is,  a  ramp  angle  of  6  degrees  and  an  initial  velocity  of  160  feet  per  second, 

A  comparison  of  Figs.  8  and  19  shows  that  in  the  C-16  crash,  maximum  ac¬ 
celerations  and  time  after  impact  at  which  they  occurred  were  somewhat  less 
than  in  the  Constellation  crash.  The  general  trend  of  the  acceleration-time 
histories  was,  however,  much  the  same. 

4.3.3  Ramp  Angle 

Fig.  20  gives  the  variation  with  ramp  angle  of  peak  longitudinal  ac¬ 
celeration  at  the  airplane  center  of  gravity  and  time  at  which  it  occurs.  This  « 

peak  is  not  the  maximum  peak  G’s  at  an  approximate  time  after  impact  of  0.  22 
seconds,  but  rather  is  the  peak  caused  by  the  initial  impact  force  at  the  nose. 

Although  results  for  ramp  angles  in  excess  of  ten  degrees  are  not  shown,  it  is 

expected  that  peak  values  would  rise  sharply  as  the  ramp  angle  increases.  • 

Figures  21  and  22  show  the  acceleration  time  histories  for  the  normal  and 

longitudinal  accelerations  for  three  ramp  angle  conditions. 

The  position  of  the  airplane  relative  to  the  base  of  the  ramp  at  impact 
is  also  important  in  determining  peak  accelerations  and  the  time  at  which  they 
occur.  For  example,  if  the  airplane  impacts  the  ramp  in  Position  B,  Fig.  21, 
as  is  the  case  in  the  Lockheed  1649  crash,  the  acceleration  time  histories  will 
be  substantially  different  from  those  in  which  the  airplane  impacts  the  ramp  in 
Position  A,  Fig.  23.  In  Position  A  the  airplane  must  pitch  through  an  angle 
equal  to  that  of  the  ramp  before  the  under  structure  of  the  fuselage  contacts  the 
ramp’s  surface,  whereas  in  Position  B,  only  a  small  pitch  angle  is  required  be¬ 
fore  the  tail  begins  to  pick  up  load.  The  fuselage  above  the  cabin  floor  will  be 
more  likely  to  fail  in  the  impact  condition  of  Position  B  since  the  fuselage  bend¬ 
ing  moment  is  considerably  greater. 

4.3.4  Friction 

The  total  frietion  force  as  given  by  Equation  (8)  of  Section  2. 1.  2  is  a 
function  of  the  normal  force  on  the  airplane  and  of  the  friction  coefficient. 

Maximum  friction  forces  occur,  then,  at  the  same  time  as  do  maximum  nor¬ 
mal  forces.  Since  peak  normal  forces  at  the  airplane  center  of  gravity  may  be 
of  the  order  of  10  to  20  times  the  weight  of  the  aircraft,  the  friction  force  is 
a  primary  quantity  in  the  determination  of  peak  longitudinal  accelerations. 

For  this  analysis,  a  nominal  friction  coefficienl  of  0.  3  was  assumed. 

This  is  the  friction  coefficient  of  aluminum  sliding  on  clay  as  given  by  Ref.  9, 
and  is  believed  to  be  close  to  that  which  will  be  applicable  to  the  Lockheed 

crash.  The  friction  coefficient  was  varied  in  the  range  of  0.  3  to  0.5,  and  the  • 

results  are  shown  in  Fig.  24.  As  expected,  the  peak  longitudinal  accelerations 
increased  with  increasing  friction  coefficient,  but  the  time  to  peak  did  not 
change. 

4.3.5  Plowing  Force 

The  plowing  force  is  analytically  determined  by  Equation  (9)  of  Section 
2.  1.  2  and  is  a  function  of  the  plowing  coefficient  and  of  the  airplane  plowing 
area.  Both  o,  these  quantities  are  difficult  to  determine,  and  particularly  the 
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IMPACT  VELOCITY  •  180  H7SEC 


FIG.  20  EFFECT  OF  RAMP  ANGLE  ON  PEAK  ACCELERATIONS 


RAMP  ANGLE  -  DEGREES 


FIG.  21  RAMP  ANGLE  VARIATION  NORMAL  ACCELERATION 
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FIG.  23  POSITION  OF  THE  AIRPLANE  AT  IMPACT 
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FRICTION  COEFFICIENT 


i  it  ltd  since  it  is  dependent  primarily  on  the  depth  of  penetration  of  the  airplane 
into  the  ground.  This,  in  turn,  is  a  function  of  the  weight  and  geometry  of  the 
airplane,  the  deformation  characteristics  of  both  the  airplane  and  ground,  and 
the  angle  at  which  the  aircraft  impacts.  However,  results  have  indicated  that 
in  the  1649  crash  peak  longitudinal  accelerations ,  which  are  most  affected  Dy 
the  plowing  force  in  low  angle  of  impact  crashes,  are  determined  primarily  by 
the  friction  force.  The  effects  of  an  increase  in  plowing  force  are  illustrated 
in  Fig.  25.  Tripling  the  plowing  coefficient  increased  the  peak  longitudinal  ac¬ 
celerations  from  2.7  G’s  to  3.4  G’s,  a  substantially  lower  rate  of  increase 
than  would  be  obtained  by  tripling  the  friction  coefficient  as  indicated  by  Fig. 

22. 


No  experimental  data  on  the  plowing  coefficient  at  the  crash  test  site 
was  available,  and  a  nominal  value  of  80  pounds  per  square  inch  of  plowing  area 
was  chosen.  This  is  similar  to  that  used  in  the  calculations  of  Ref.  9.  A  nom¬ 
inal  plowing  area  of  700  square  inches  was  estimated  based  on  the  data  of  Ref. 

9  and  a  study  of  the  1649  geometry. 

4.3.6  Damping  Forces 

Damping  forces  are  considered  to  be  those  forces  in  the  crash  which 
are  directly  proportional  to  the  velocity  of  the  airplane  normal  and  downward 
into  the  ramp.  They  arise  primarily  from  the  friction  of  deformation  of  the 
ground  and  uneerstructure  of  the  airplane.  For  the  1649  crash  a  structural 
damping  of  0.25  was  assumed. 

The  damping  was  found  to  have  seme  effect  on  the  acceleration  time 
histories.  For  this  reason,  and  since  the  structural  damping  coefficient  is 
impossible  to  determine  with  any  great  amount  of  precision,  it  was  varied  in 
the  range  of  0. 100  to  0.500.  A  summary  plot  showing  peak  normal  accelera¬ 
tions  at  the  center  of  gravity  for  each  of  three  damping  coefficients  is  given  in 
Fig.  26. 


In  general,  a  decrease  in  damping  coefficient  was  found  to  slightly  in¬ 
crease  the  peak  normal  accelerations  at  0.  07  seconds  after  impact  and  decrease 
both  the  magnitude  and  time  of  the  maximum  peak  accelerations  occurring  at 
approximately  0.24  seconds.  This  correspondingly  reduced  maximum  longi¬ 
tudinal  accelerations. 

4.3.7  Ground  Elasticity 

Inclusion  of  the  effective  spring  stiffness  of  the  ground  had  little  effect 
on  the  results.  The  ground  spring  was  considered  to  act  in  series  with  the 
structural  spring  of  the  understructure  of  the  fuselage  as  per  Section  2.  1.  2. 
Experimental  data  from  the  crash  test  site  gave  a  California  Bearing  Ratio 
(CBR)  rating  of  30-50  to  the  soil  in  the  impact  area.  This  was  converted  to  an 
approximate  spring  stiffness  of  500,000  pounds  per  inch  which  was  substanti¬ 
ally  greater  than  that  of  the  structural  springs.  For  this  reason,  when  the 
ground  spring  was  placed  in  series  with  the  structural  spring,  the  stiffness  of 
the  structural  spring  dominated  the  net  resultant  stiffness,  and  the  effect  of  the 
ground  spring  was  negligible. 
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IMPACT  VELOCITY  -  180  FT/SEC 


STRUCTURAL  DAMPING  COEFFICIENT 


5. 0  CONCLUSIONS 


•  At  initial  impact  of  the  Lockheed  Model  1649  at  a  velocity  of  180  feet  per 
second,  a  force  from  the  ground  acting  at  the  nose  of  the  airplane  'ill 
cause  it  to  pitch  up  sharply, 

•  During  this  initial  impact  phase,  peak  normal  accelerations  of  11,  0,  and 
-3  Gs  will  be  developed  approximately  0.  03  seconds  after  impact  at  Body 
Stations  180,  682,  and  1176  respectively.  At  this  point  in  time,  an  air¬ 
plane  longitudinal  acceleration  of  1  G’s  will  be  experiencr  d. 

•  The  nose  of  the  airplane  will  bend  upwards  10  inches  relative  to  the  center 
section  of  the  fuselage  at  a  time  after  impact  of  0. 14  seconds.  This  de¬ 
flection  will  probably  be  of  sufficient  magnitude  to  cause  the  ultimate 
strength  of  the  fuselage  above  the  cabin  floor  to  be  exceeded. 

•  If  the  fuselage  above  the  cabin  floor  should  fail,  all  analytical  results  be¬ 
yond  the  time  of  failure  will  be  questionable  to  a  degree  dependent  on  the 
type  of  failure  that  occurs. 

•  If  the  fuselage  above  the  cabin  floor  does  not  fail,  maximum  normal  and 
longitudinal  accelerations  during  the  6  degree  ramp  crash  will  occur  at 
a  time  after  impact  of  0.24  seconds.  Maximum  normal  acceleration  at 
Body  Stations  180,  682,  and  1176  will  be  -17,  8,  and  35  G’s  respectively. 
Maximum  longitudinal  accelerations  will  be  4  G’s. 
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6. 0  RECOMMENDATIONS 


Further  studies  are  necessary  to  better  define  the  effects  of  initial  crash 
conditions,  and  airplane  weight,  geometry,  and  structural  characteristics 
on  crash  loads  for  different  aircraft. 

Improved  analytical  or  experimental  data  are  required  for  determining 
the  crushing  strength,  plowing  area,  and  structural  damping  of  an  airplane 
during  a  crash. 

The  influence  of  airplane  roll,  yaw,  and  lateral  translation  on  crash  loads 
should  be  studied. 

If  a  sufficient  amount  of  crash  loads  data  can  be  obtained,  it  should  be 
used  to  improve  the  accuracy  of  required  airplane  crash  load  factors.  It 
should  also  be  used  to  define  dynamic,  rather  than  static,  test  specifi¬ 
cations  for  items  of  equipment  which  must  function  during  a  crash.  In¬ 
cluded  in  this  equipment  would  be  passenger,  crew,  and  cargo  restraints, 
and  shock  absorption  systems. 
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